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NMSSM scenarios are investigated to explain an excess in the opposite-sign dilepton mass distri¬ 
bution in events with dilepton, jets and missing transverse energy reported by the CMS experiment. 
We show that the NMSSM scenarios can possess unique features to explain this excess, and can 
be distinguished from the MSSM scenarios in the ongoing LHC runs as well as direct detection 
experiments. 


I. INTRODUCTION 

Recently, CMS reported of an excess of lepton pairs [T] 
with energy below the Z mass, in a final state of 
l~^l~jj + pT where ^ denotes the missing transverse 
energy. Beyond the the Standard Model (BSM) theoret¬ 
ical attempts have been made to account for such an ex¬ 
cess in the context of the Minimal Supersymmetric Stan¬ 
dard Model (MSSM) [IH3] using cascade decays from the 
sbottom pair production at the LHC, a leptoquark sce¬ 
nario |3] and superstring inspired models [3, where the 
cascade decay of new particle states give rise to a pair of 
leptons, two associated jets, and new invisible particles 
that escape the detector as pr- 

Two same-flavor, opposite-sign leptons can be pro¬ 
duced from a cascade decay that has particles decay 
into leptons in the intermediate state during the process. 
Such lepton partners and heavier states, which give rise 
to the cascades, are readily available in supersymmetric 
models. The sbottom pair production and its cascade 
decays into the next to lightest neutralino which sub¬ 
sequently decays into two leptons and the lightest neu¬ 
tralino via an intermediate slepton state is a very inter¬ 
esting option to explain the excess. Two b jets are also 
produced along with l~^l~ in the final states. The 

existence of two b jets in the signal provides an inter¬ 
esting prediction arising from this scenario which will be 
checked in Run II. Prior to this new result, l~^l' jj + pT 
was considered as a possible final state from the stop 
decay [B]. Since we have i —?> i + X 2 ) expect lep- 
ton(s) from top decay in addition to X 2 This 

doesnt support the GMS edge paper. In this paper we 
focus on a well motivated Next-to Minimal Supersym¬ 
metric Standard Model [2] (NMSSM), which introduces 
an additional singlet superfield into the MSSM. The ob¬ 
served Higgs mass at the LHG [5] can be accommodated 
naturally if the coupling A between the singlet and the su¬ 
persymmetric Higgs fields is large and this new term also 
provides a solution to the /i problem of the MSSM (see [5] 
and references therein). Recently, signals of such a sce¬ 
nario at the LHC were investigated and possible ways to 
distinguish from the MSSM were also discussed [TO] . 



FIG. 1: Pair production of sbottoms and one of its cascade 
decay chains can give rise to a signal of lepton pairs with en¬ 
ergy below the neutralino mass difference. The other sbottom 
in the diagram presumably decays directly to 

The singlet superfield in the NMSSM gives rise to a 
new neutralino (singlino) to the gaugino sector, besides 
other modification on the MSSM particle spectrum. This 
paper investigates the effect of the singlino that creates 
more freedom to realize the aforementioned cascade de¬ 
cays of sbottom, which alleviates a relatively tight re¬ 
quirement on the MSSM gaugino mixings. Further, since 
the neutralino sector is modified the NMSSM explanation 
will be associated with distinguishable predictions. The 
cascade process to explain the endpoint in the context of 
the NMSSM is discussed in Section |TT] A collider analysis 
is carried out in Section Hill that examines the NMSSM’s 
explanation of the CMS dilepton excess. We discuss our 
results and conclude in Section m 

II. SQUARK-NEUTRALINO CASCADES 

The CMS excess in the relatively low energy ‘endpoint’ 
in opposite sign dileptons can rise from the cascade de¬ 
cays of one or more next-to-lightest supersymmetric par¬ 
ticles (NLSPs) in the supersymmetric neutralino sector, 
where the second lightest neutralino can decay via a two- 
step process, 

X 2 ^ ^ Xli ( 1 ) 

that yields a pair of leptons of the same flavor, whose to¬ 
tal energy distribution is limited by the mass difference 
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M^o — M^o {I denotes both both left and right sleptons of 
the first two lepton flavors.). Thus a neutralino spectrum 
with M^o ~ M^o +70 GeV can lead to such a dilepton ex¬ 
cess if the slepton masses lie between the two neutralino 
masses and if the second neutralino can be produced at 
adequate rates. 

The QCD-dominated production of sbottom^ can lead 
to the correct pair-production rates for the CMS exper¬ 
iment. A mostly right-handed sbottom (denoted by bi 
here) can be a perfect candidate [T] that decays into the 
required neutralinos. The Feynman diagram for the pro¬ 
duction process is illustrated in Fig. The sbottom 
is preferably right-handed as it does not require a left- 
handed stop to be equally light in mass. Moreover, the 
CMS non-observation of more than two leptons in this 
channel puts a constraint on the decay of such a sbottom, 
i.e., the branching fraction, BF( 6 i —)• Xib), should domi¬ 
nate over BF( 6 i —>• x^b) unless the latter produces many 
lepton pairs in a final state from the pair of sbottoms. 
This requirement, however, may force certain relations 
among the neutralino mixings. 

In the MSSM, a sbottom can decay into a b quark and 
a neutralino via either the U(l)y gauge coupling or the 
Yukawa coupling. The U(l)v gauge coupling depends 
on the Bino (B) component of the daughter neutralino, 
while the Yukawa coupling depends on the down-type 
Higgsino (Hd) component and the size of the coupling yd- 
Since the MSSM’s LSP needs to be dominantly Bino to 
avoid increasingly severe experimental constraints (from 
indirect mnn] and direct [13] detections). However, the 
Bino content of the LSP would dictate the sbottom decay 
(by a few orders of magnitude) and leaves with small B 
for X 2 for the hi x^h decay when M^o — M^o « 70 GeV. 

The BF to bi —>■ x^^ can be increased when tan /3 is large 
which causes larger Yukawa coupling contribution [2|. 

In the NMSSM, however, the option of a singlino LSP 
opens up alternative gaugino mixing scenarios to realize 
sbottom decays into both Xi and X 2 ia a suitable way, 
which differ from the MSSM. Consequently, the collider 
signals arising from the direct production of gauginos will 
be helpful to discern MSSM from NMSSM, 

(A) We can satisfy the CMS excess for mostly singlino 
LSP and a heavy Bino. The components of B in both 
AijA 2 are small but comparable in magnitude for heavier 
B. This allows a small tan/3 scheme where both bi —>■ 
Xib,X 2 ^ decays occur via the U(l)y coupling, while the 
down-type Yukawa contribution to the decays is small. 

The NLSP can be wino and/or Higgsino. A wino- 
NLSP case is shown in Table |T| In comparison, when the 
NLSP is mainly Higgsino, X 3 would be relatively light 


^ The dilepton mass endpoint can also be produced from the decay 
of a light stop, but when stop is pair produced we also expect 
to see multilepton final state and the non-observation of more 
than two leptons along with 2 jet and pi makes this scenario 
not preferred. 


and the sbottom may decay into two NLSPs. However, 
the Higgsino-NLSP case tends to allow the LSP to have 
a larger Higgsino mixing which faces constraint from the 
direct detection result unless we chose the correct sign of 
gaugino, Higgsino mass parameters to cancel the Higgs 
contribution in the direct detection amplitude M- 
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TABLE I: The NMSSM benchmark points that yield same- 
flavor opposite-sign lepton pairs in sbottom cascade decays. 
The mass spectrum is evaluated using NMSSMTools [15] and 
its values are given in GeV. I denotes both both left and right 
sleptons of the first two lepton flavors. 

(B) We can satisfy the CMS excess for mostly Bino 
type Xi, where X2 mostly singlino and Xs 4 consist 
of mostly Higgsinos and Wino. Just like the MSSM, a 
large tan /3 is required to boost the decay via down-type 
Yukawa coupling to the Hd component in X 2 ^ which is 
closer in mass to the Higgsinos in comparison to the much 
lighter x^. The major difference of this scenario from the 
MSSM is that the lightest chargino xt mass is close to 
X 3 rather than X 2 ^ consequently can be heavier than 
that of the MSSM. This allows a wider mass range of the 
sbottom after satisfying — M-+ < Mt so that bi does 
not decay into top quarks. 

(C) We can also satisfy the CMS excess for a mostly 

singlino type Xi where X 2 i® mostly Bino. Since X 2 i® 
mostly Bino, BF(&i —>■ X 2 ^) i® large and yields a large 
number of final state leptons via sleptons situated in be¬ 
tween the two neutralinos which may not be a suitable 
option. A virtual slepton mediated three-body decay 
{X 2 1*1^,I* —f however, can give the correct 

21 + pT rate. This parameter space is represented by 
point C which shows that the slepton masses are much 
higher than the Xi 2 mass range, e.g., the slepton masses 
are almost at TeV scale for point C. It is interesting to 
note that even if we imagine a scenario where the slepton 
masses are very close to either X 2 or Xii ^-nd only one of 
the leptons from each sbottom cascade is visible, the in¬ 
variant mass of the two leptons from different cascades 
can easily be more than 70 GeV and contradicts with 
observation. 


III. COLLIDER SIGNAL 

Here we discuss the dilepton yield from the NMSSM 
benchmark points at the 8 TeV LHC. At each point, a 
mass gap is kept at 70 GeV between the two lightest 
neutralinos that limits the energy of the lepton pair. X 3 is 
above bi and all sfermions other than I and bi have multi- 
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TeV masses and decouple from our study^. The lightest 
chargino xt mass is not lighter than —Mt, which can 
be used as one aspect that the NMSSM spectrum differs 
from the MSSM’s. The bi is mostly right-handed and its 
decay BFs are listed in Table [ill together with the BF of 
X 2 ’s decay to I and the lepton reconstruction efficiency. 
The BF([ —>■ Ixi) decay is 100% for all of our benchmark 
points, since in all these cases the LSP is either a singlino 
or bino. 







A 

69% 

31% 

25% 

38% 

B 

90% 

10% 

57% 

38% 

C 

10% 

89% 

11% 

26% 


TABLE II: hi decay branching fractions and the dilepton 
selection efficiencies for the benchmark points, eu is given 
in Eq. for point C in column 4 should read as 


We generate inclusive sbottom pair production events 
with 0-2 associated jets at 8 TeV in MadGraphS 
v2.2.2 [T^ with CTEQ6.6 [T7j parton distribution func¬ 
tions. Pythia v6.426 [15] is used for showering and 
hadronization and PGS4 m is used for detector sim¬ 
ulation in which the electron and muon detection effi¬ 
ciencies are assumed to be 92% and 98%, respectively. 
To avoid double-counting of jets the MLM jet matching 
scheme |5D| is implemented. 

The LO production cross section ctlo ~ 500 fb for 
sbottom with approximately 360 GeV mass. This is 
scaled up by a K-factor of 1.7 to the NLO value that is 
obtained from the package Prospino [21] . The dilepton 
signal rate is then. 


(20 fb • ctnlo • Aeff, (2) 


where Aeff denotes the total event selection acceptance. 
The selection criteria for the jet and lepton objects are as 
follows: PT(jet)> 40 GeV with \r]\ < 3, p 7 ’(lepton)> 20 
GeV with \r]\ < 2.4, excluding 1.4 < \ri\ < 1.6; the central 
region is defined as \r]\ < 1.4. 


^ We remain agnostic about the exact mass of other heavy squarks 
and ignore the small change in the NLO b production cross- 
section due to variation of them. 


Event selection 

Relative Efficiency 


A 

B 

C 

cnlo (fb) 

660 

854 

854 

Nj > 2(3) + pr> 150(100) GeV 

32% 

37% 

25% 

Two isolated OSSF leptons 

3.2% 

2.3% 

3.2% 

Dileptons in the central region 

85% 

85% 

84% 

Overall acceptance 

0.85% 

0.72% 

0.66% 

Number of events at 20 fb~^ 

112 

124 

112 


TABLE III: The selection efficiencies for the I'^l jj + f/r sig¬ 
nal at benchmark points A, B and C. 

The selection efficiency flows for points A, B and G 
are shown in Table m Overall acceptances are similar 
for A and B. The overall acceptance for point G is small 
due to the fact. The first cut is the “OR” cut from the 
CMS analysis [Tj: it selects events with at least 2 jets 
and pT >150 GeV, or at least 3 jets and pr >100 GeV. 
We find that point C shows lower efficiency at this stage 
compared to points A and B. This is due to the fact that 
bi has a larger branching ratio into bi —> b+pl for points 
A and B where i® mostly Bino compared to the point 
C where xj i® mostly singlino. A large mass gap 
120 GeV) between p and Xi causes the existence more 
higher px jets for points A and B. On the other hand, 
for point C, bi has larger branching ratio into b + X 2 
and a smaller mass gap (~ 50 GeV) between bi and X 2 
produces lower px jet for the point C which causes lower 
efficiency after we require pT(jet)>40 GeV. The fraction 
of isolated leptons is a combination of the 5, X 2 decay 
branching fractions and dilepton selection efficiency eu, 

2-x{l-x)eu, 

x = mbi^X°2b)BF{x°2^ll). 

eu is a fraction of events with at least two reconstructed 
OSSF leptons passing the lepton selection criteria out 
of the events with two OSSF leptons from the X 2 de¬ 
cays. In principle, eu can vary with the lepton energy 
that is determined by the mass difference between the 
sleptons and neutralinos. At points A and B, when 
Mj ~ + M^o)/2 with both leptons at Ei = 35 GeV, 

we get maximal acceptance of eu « 38%. At point G 
the leptons arise from a three body decay and the energy 
partition becomes uneven, and the acceptance can suffer 
if one of the leptons is too soft. The lepton pairs are 
further required to be in the central region with pseu¬ 
dorapidity \r]\ < 1.4, and finally the overall acceptance 
gives the accumulated selection efficiency. 

At all three benchmark points, a dilepton signal of 
around 120 events are obtained. The dilepton mass 
distributions for points A and B are similar, and they 
both show a clear dilepton invariant mass endpoint at 
M^o — M^o as each of two leptons has a fixed 35 GeV 
energy, where the invariant mass of the dilepton maxi¬ 
mizes. The endpoint for point C is less pronounced due 
to the possible unequal energy between the two leptons 
in the three body decay of X 2 - In Fig-[^ we show the dis- 
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FIG. 2: Invariant mass distribution of the opposite sign same 
flavor dileptons after implementing all selection cuts. The dis¬ 
tributions are based on MC sample of 1000 fb~^ integrated 
luminosity but normalized to 20 and 100 fb“^ integrated lu¬ 
minosities. 

tributions for points A and C. The distribution for point 
B is identical to that of point A and is not shown. As the 
lepton energy is well measured at the LHC, it is possible 
to discriminate point C from other points by performing 
shape analysis with finer binning and higher integrated 
luminosity, as is shown in the inset. The statistical fluc¬ 
tuations in Fig. [^corresponds to a MC sample equivalent 
to 1000 fb“^ integrated luminosity. All three benchmark 
points have less than 4% multilepton final states com¬ 
pared to the dilepton final states as shown in Table |TV| 
We have chosen px > 10 GeV for the third and fourth 
lepton in the multilepton final states. 


Number of multilepton events at 20 fb“^ 

A 

B 

C 

Three lepton 

2.14 

2.10 

4.46 

Four leptons 

1.10 

0.7 

1.52 


TABLE IV: Number of and Aljj + fk events at benchmark 
points A, B and C. 


IV. DISCUSSION AND CONCLUSION 

The dilepton excess, if confirmed at the ongoing run of 
the LHC, will be able to probe SUSY models by investi¬ 
gating other signals associated with the models. In the 
case of the MSSM, the dilepton excess requires the exis¬ 
tence of a Higgsino-type second-lightest neutralino and a 
chargino for masses around 300 GeV. 14 TeV LHC can 
produce these Higgsinos directly. However, in such a sce¬ 
nario, the direct dark matter detection cross-section can 
potentially be a problem due to an appreciable mixture 
of Bino and Higgsino in the LSP which would push us to 
use the negative sign of p. This problem can be avoided 
in the NMSSM which also exhibits other distinguishable 
features which we discuss below. 

For point A, the NLSP can be mostly wino. In this 
case, it is possible to produce winos directly and the 


production cross-sections by Drell-Yan and vector bo¬ 
son fusion processes are larger by almost a factor of 4 
compared to that in the Higgsino NLSP scenarios. The 
current bound(~ 300 GeV [SSIEl]) still allows this point 
since the mass difference between the , X 2 Xi 
small and X 2 H decay branching is less than what is 
assumed by the CMS and ATLAS. The ongoing run will 
probe this point which does not suffer from the direct 
detection constraint. 

For point B, more neutralino states are involved to 
explain the excess. In this case the Higgsinos, which 
constitute mostly the Xg 4 , are lighter than winos and 
they can be also be produced directly at the LHC. The 
350 GeV Xg {m^o is 410 GeV) has a 75% decay BF into 
Xg —t .^Xi most of the remaining BF is into rr. The 

Xi is at 333 GeV and has appreciable a BF into lU + Xi- 
We have already shown that the sbottom decays mostly 
into X 2 subsequently into x°, whereas Xa 4 and xt 
decays predominantly to Xi- The direct productions of 
X 2 and Xi are small since they are mostly singlino and 
Bino respectively. Therefore, at least four neutralinos 
and one chargino will be potentially discovered at the 
ongoing run from two different production processes and 
can be used in a complimentary way to check whether 
this model point explanation for the excess is correct. In 
addition, just like the Bino-Higgsino case in the MSSM, 
this scenario suffers from the direct detection constraints. 
However, flipping the sign of the gaugino and Higgsino 
mass parameters, can suppress the nucleon-LSP scatter¬ 
ing. 

For both points A, B and the MSSM, the mass-gap 
needed for these points to explain the excess is about 35 
GeV. It can be difficult to probe the sleptons directly in 
the current runs of the LHG due to relatively small mass- 
gap between I and the LSP. Even boosting the system 
with additional jets does not help, since such searches lose 
their efficacy for a mass splitting above ~ 25 GeV [25] . 

In contrast, at point G, the sleptons are relatively 
heavy in order not to overproduce event with more than 
two leptons. If the sleptons are found to be a lot heav¬ 
ier than the NLSP, then the point G will be needed to 
explain the GMS excess. However, the direct slepton pro¬ 
duction can be probed up to ~ 700 GeV [121 [IS] at 3000 
fb“^ of integrated luminosity. However, since X 2 Xi 
are mostly Bino and singlino respectively, the direct pro¬ 
duction cross-sections of them are very small. Further, 
unlike points A, B and MSSM, this point does not have 
any light chargino associated with the lightest neutrali¬ 
nos which can be probed at the LHG. Therefore, if no 
other lighter neutralinos and charginos are found but the 
endpoint in the dilepton distribution still persists, the 
point G will provide the explanation. This point also 
does not suffer from the direct detection problem. It is 
also possible to discriminate point G from other points 
by performing shape analysis. 

In conclusion, if the GMS excess is proved to be cor¬ 
rect at the ongoing run of the LHG, it will be possible 
to find an explanation in the context of NMSSM, which 
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is not ruled out by the direct detection experiments and 
with unique features in the neutralino sector compared to 
the MSSM. It will be feasible to establish these NMSSM 
scenarios by investigating the direct productions of neu- 
tralinos, charginos and sleptons. 
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